Orbital Magnetism in Ensembles of Parabolic Potentials 
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We study the magnetic susceptibility of an ensemble of non-interacting electrons confined by 
parabolic potentials and subjected to a perpendicular magnetic field at finite temperatures. We 
show that the behavior of the average susceptibility is qualitatively different from that of billiards. 
When averaged over the Fermi energy the susceptibility exhibits a large paramagnetic response only 
at certain special field values, corresponding to comensurate classical frequencies, being negligible 
elsewhere. We derive approximate analytical formulae for the susceptibility and compare the results 
with numerical calculations. 
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The interest in the magnetic properties of ensembles 
of mesoscopic systems has increased considerably in the 
last years • The main motivation for the theoretical 
investigations recently carried out is the experimental re- 
sults obtained by Levy et al 0] for an ensemble of square 
billiards in the ballistic regime. It is now understood 
that, when averaged over a large ensemble of similar sys- 
tems, the magnetic susceptibility x of regular systems is 
enhanced with respect to the Landau susceptibility xl 
due to the coherent contribution of families of periodic 
orbits. For chaotic systems x is usually small, of the or- 
der of xl, but bifurcations might play an important role 
in increasing the susceptibility ||. Also, for square bil- 
liards, the averaged susceptibility is always paramagnetic 
at low magnetic fields. This behavior seems to be generic 
of regular billiards . 

The purpose of this letter is to study the magnetic 
susceptibility of an ensemble of noninteracting two- 
dimensional electron gas confined by parabolic quantum 
wells at finite temperatures. We show that the behav- 
ior of x as a function of the magnetic field B and Fermi 
energy p is very different from that of billiards when en- 
semble averages are considered. The main reason for the 
this strong difference resides on the resonances exhibited 
by the system as the magnetic field is varied. The sus- 
ceptibility of a single harmonic oscillator in a magnetic 
field has been considered before by Prado et al H and 
Nemeth [Q. In this work we derive simple analytic ex- 
pressions for x which are directly amenable of ensemble 
averages. 

In the canonical ensemble the magnetic susceptibility 
per particle x — ~ (1M0 d 2 F/dB 2 measures the sensi- 
tivity of the Helmholtz free energy F to the magnetic 
field B. For a non- interacting model F can be computed 
exactly if the number of particles N and the temperature 
T are not too large ||. However, for large N and temper- 
atures of the order of the the mean level spacing (in units 
of Boltzmann constant fes), X can be computed quite ac- 
curately by using an auxiliary grand-canonical potential 
where the average number of particles is kept fixed by 
properly adjusting the chemical potential for each value 
of the magnetic field M . The grand-canonical potential 



is given by 



V=-J 3 J dEp(E) ln(l 



„P(p-E) 



(1) 



where p is the density of states, (3 = 1/fcsT and p is 
the chemical potential. In the semiclassical limit p can 
be separated into a mean term p° plus oscillatory con- 
tributions p osc that are usually written in terms of pe- 
riod orbits 10 1. After substituting p = p a + p osc in 
Eq. ([j]) , the grand-canonical potential also separates into 
V°(fi) + V osc (p). Following Ullmo et al @ we define a 
mean chemical potential p° from N — J dE p(E)f(E — 
p) = J dE p°(E)f(E - p°) where f(x) = 1/(1 + e^) is 
the Fermi-Dirac distribution function. Using the thermo- 
dynamical relation F(N) = V(p) + pN with p = p(N) 
obtained from the equations above, it can be shown j^] 
that, in the semiclassical limit, the Free Energy can be 
written as a sum of three terms, F = F° + AF 1 + AF 2 , 
where F° = V°(p°) + p°N does not depend on the mag- 
netic field B, AF 1 = V osc (p°) and 



2p°(p°) 



dE p osc (E)f(E - p°) 



-I 2 
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In the case of a parabolic confinement, the Hamiltonian 
of an electron of charge e and effective mass m* can 
be written directly in terms of action variables as |ll]] 
H(I) = n ■ I with (\ = + «2 = - £-) 
and £± = \J (u)\ ± lo-i) 2 + e 2 B 2 /m* 2 . The quantum me- 
chanical energy levels are therefore given by Ek 1 k 2 — 
hfli(ki + 1/2) + hVt 2 {k 2 + 1/2). To compute the p° and 
posc we f u ow Berry and Tabor jl2j and write the semi- 
classical density of states as 



1/2)) 



(3) 



where m = (mi, 7712) and the factor 2 takes care of 
spin degeneracy. Using the Poisson sum formula we find 
p°{E) = 2E/(h 2 uj 1 uj 2 ) and 



-iTr(m 1 +m 2 ) 
e Pn 



(4) 
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with 



p m = I d£ 



3 27rim-I(£)/7» 



(5) 



where £ varies along the energy surface in the (/i,^) 
plane of the action variables. The prime above the sum- 
mation symbol means that the term (mi, 7712) = (0,0) 
(responsible for p°) is excluded. Explicitly we get I\ — 
E/VLx - ffi 2 /Q and I 2 = CliZ/tt for < £ < M2/(fiifi 2 ) 
with n = y/ttl + Of. For generic Hamiltonians the inte- 
gration over £ can performed within the stationary phase 
approximation and gives a semiclassical expression for 
p osc . In the present case the integration can be carried 
out exactly and results in 



' 2 (_i r i+™2 e Hffrni< m 2"i+"i"2) 

A^m 7rh(m2^i —mi Q2) 



xsin [?I^( TO 2^i -mxQ 2 )} 



(6) 



Substituting the above expression for p osc in the formulas 
for AF 1 and AF 2 we get 



af 1 = y 



(-1)' 



S, 2 wiw 2 (72-7i) /3 



sin (7 2 )(l-(5 m2 o) sin (7i)(l-(5 mi0 ) 



72 sinh 



(7) 



and 



AF 2 - J- iV" „(- 1 ) mi 



(72-71) P 



COS (72)(l~l?m 2 0) _ COS (7l)(l-^ mi o) 

sinh ( 



sinh (^2.) 



(8) 



The double prime in the summations means that only 
the integers (7711,777,2) in the upper half plane, minus the 
negative m\ axis, are included. We have also defined 

7i = 27T/J°77l;/(?iSlj). 

Both expressions for AF 1 and AF 2 can be simplified 
if one notes that, as the magnetic field is varied, the ratio 
Oi/ O2 passes densely through rational numbers, where 
all the classical orbits of the system are periodic. To those 
field values there also corresponds a large degeneracy of 
the energy levels. We therefore restrict our attention ini- 
tially to the neighborhood of these values of B only. Let 
B = B nm be such that f^i = 77SI0 and fi 2 = 777^0 , i.e, 
^1/^2 = 77/777 and fio = yj "* t ^ 2 is the frequency of the 
classical periodic orbits. At those points the denomina- 
tor in (^) and (JsJ) vanishes for all (7771,7772) of the form 
(pn,pm) for all p > and we see that the main contribu- 
tions to AF 1 and AF 2 come from these resonant terms. 
Defining the function 



S(B) = 



sin [C(B)} 
C{B) 



(9) 
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FIG. 1. The first 290 single particle energy levels normal- 
ized to e = E/(hui) as a function of the magnetic field B. 
The triangles on top indicate the resonances with m = 1 and 
n = 2, ...,20. 



where C(B) = ±( 72 
at B nm , 



mtti) vanishes 



aw* 7i J = 1^7(^2 ■ 
t , we get, after rearranging the trigonometric func- 
tions and considering only the term p — 1 , the following 
approximated expressions: 



AF 1 



(-l)"+-/i°F(/3) ,27^° qm , 

7T Z 77?77 niln 



and 



AF 



, sin 2 ( 2 ^°^ 



) S 2 (B) 



(10) 
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where R(/3) = (2tt 2 ')/ '(fifffio) sinh -1 [(2tt 2 )/(^/3O )] is 
a temperature dependent factor that diminishes expo- 
nentially the susceptibility for large T's and large Oo's. 
The dependence of the free energy on the magnetic field 
has been reduced to S(B) and the magnetic suscep- 
tibility can be readily computed. We recall that the 
Landau susceptibility for the oscillator is given by 0] 

XL = -(e/777*) 2 ?7 2 /(6 M ). 

The expressions just derived describe the behavior of 
a single system in the thermodynamic limit. In the ex- 
periment with square billiards of Levy et al p|, however, 
only the average properties of an ensemble of systems 
were measured. The individual members of the ensemble, 
although very similar, present small differences among 
themselves. Besides, the number of particles confined in 
each of them might vary slightly. To account for these 
fluctuations further averages have to be performed . 
As in the case of billiards ]jj the oscillatory contribution 
of AF 1 to x vanishes under an average over the Fermi 
energy, or number of particles, for dispersions Sp, of the 
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FIG. 2. Magnetic susceptibility % as a function of /x° at 
73 = and T = 0.3 mK. The inset shows x as a function 
of B for a single system, without any average, for the same 
temperature and TV = 500. In both cases the full line rep- 
resents the numerical exact result and the dotted line shows 
the analytical calculation (\ L + X 1 + X 2 )- 
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FIG. 3. Susceptibility averaged over N as a function of the 
magnetic field for three different temperatures. In all cases 
N = 500, SN/N = 0.2, the full line represents the exact 
result, the dotted line shows the resonant approximation and 
the dashed line shows the full analytic formula, Eq.i 



order of HQo ■ The contribution of AF 2 remains for 
the parabolic potential as it does in the case of billiards. 
The main difference here is that, due to the fact that the 
density of states increases linearly with the energy, the 
contribution of AF 2 is of the same order in /jP than that 
of AF 1 . Also, in terms of particle number, the relative 
dispersion SN/N necessary to kill AF 1 falls as 1/yN. 
Therefore, for large iV's, even very small dispersions will 
effectively wash out AF 1 . The resulting susceptibility, 
after performing the average, is 



(X) 



N 



2ir 2 Nnm 



R 2 (13) 



2 S 2 (B) 
dB 2 



(12) 



Therefore, since d 2 S 2 (B) / dB 2 is proportional to N and 
has a negative peak at F nm , x exhibits a positive peak 
at each resonance whose strength goes as /z°/nm. 

In what follows we present numerical calculations per- 
formed with ojx — 5.4 x 10 8 s~ 1 , w 2 = 0.9wi and m* = 
0.067m e , which is the electron effective mass for a GaAs 
quantum well. Fig.|l] shows the first 290 energy levels as 
a function of the magnetic field B. The arrows on top in- 
dicate the position of the most relevant resonances. Fig.^ 
shows x as a function of /jP at B = for a single system, 
without the average. The inset shows x as a function 
of B for 500 particles, corresponding to approximately 
fx = 7.5 x 10~ 6 eV. In both cases T = 0.3mK, the full line 
represents the numerical exact result and the dotted line 
shows the result derived from Eqsj?] and |j| The agree- 
ment between exact and analytic results is very good. 
The approximate formulas, Eqs.^ and [0] also give very 



accurate results for x close to the resonances. Fig.||shows 
(x)n, the susceptibility averaged over particle number, as 
a function of B for N = 500 and dispersion SN/N = 0.2 
computed directly from the energy levels (full line) , from 
the resonant formula ( |l2| ) (dotted line) and from the full 
analytic expression Eqs.(^,||) (dashed line) for three dif- 
ferent values of the temperature. The last two curves also 
include the contribution of xl- Once again we found an 
excellent agreement between the exact and approximate 
calculations (notice that the exact and analytical curves 
involve the separate calculation of x for various N's be- 
fore the average is performed) . 

We finally consider averages over different confining 
potentials. In order not to introduce too many param- 
eters we define each member of the ensemble to have 
lui = xujiq and LO2 — 2^20 j with a gaussian distribution 
of x around x = 1. This can be considered as a size 
average, since x changes the available area in coordinate 
space without changing the shape of the potential. In 
what follows we use an extra index to indicate quanti- 
ties computed with x = 1. Keeping the number N of par- 
ticles fixed we see that // = h\/ Nx 2 ujiquj2o = xjjP and 
£lp(B) — xQqq. Therefore, the oscillations in ( |l0| ) and 
(|TT|), which depend on the ratio (jP/Qq, are not affected 
by the average. However, since £li(B) = xQio(B/x), 
S(B) = Sq(B/x) and averaging over x is equivalent to av- 
erage over B/x. Writing x — 1 + Sx, B/x ~ B — BSx and 
we see that the average is not effective for small values of 
the magnetic field. Therefore, the resonant peaks at large 
-B's tend to be smoothed out, enhancing the susceptibility 
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FIG. 4. Size averaged susceptibility (see text) as a function 
of the magnetic field for three different temperatures. In all 
cases 8x — 0.1, the full line represents the exact result and 
the dotted and dashed lines show the contributions of the full 
analytic calculation and that of \ 2 alone respectively. 



dependency of p° with the energy, on the other hand, is 
a generic property on smooth confinements and plays an 
important role in balancing the relative contributions of 
AF 1 and AF 2 at low temperatures. We notice that the 
oscillator parameters and magnetic field can be scaled in 
order to allow for experimentally accessible values. If the 
frequencies are both multiplied by /, the Fermi energy, 
the susceptibility and the magnetic field are also multi- 
plied by /, whereas the density of states scales as 1/f. 
For / = 10 3 , for instance, we would still be considering 
fields of the order of 0.1 Tesla. 
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at the non-resonant region, close to B = for the cur- 
rent value of the parameters. The non-resonant regions 
are described approximately by Eq.(^) with (mi, m?) — 
(1,1) and show an oscillatory behavior with frequencies 
7i and 72. Expanding ji(B/x) ~ Ji{B) — B5xj t and 
imposing \ji(B/x) — Ji(B)\ — 2ir we find that the oscilla- 
tions in x die out for eB/m* ~ y/hfl(uj 2 — wf)/G"oR^I) 
where f2 is the smallest between Qi and fl2- This is 
confirmed by the numerical data displayed in Fig.^ for 
5x = 0.1 and different temperatures. Notice that the 
contribution of AF 2 can be neglected because it falls off 
much faster with the temperature than AF 1 . 

In conclusion, the magnetic response of an ensemble 
of two-dimensional electron gases confined by parabolic 
potentials is qualitatively different from that of an en- 
semble of billiards. The main features of the problem 
can be understood with the help of a resonant approx- 
imation for the density of states. For an ensemble of 
identical quantum wells, each holding slightly different 
number of electrons, the magnetic response is enhanced 
only at the resonances. When a dispersion in the size 
of the oscillators is is included, the sharp response at the 
resonances are smoothed out but the susceptibility at low 
fields remains oscillatory as a function of the Fermi en- 
ergy, and not necessarily paramagnetic as in the case of 
billiards. We emphasize that the nature of our approxi- 
mations are different from those of ref. since here it 
is validy for all values of B, not only in the limit of small 
fields. The large peaks exhibited by the susceptibility at 
the resonances, is very peculiar of the oscillator but the 
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